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ABSTRACT 

The moment of inertia of the pulsar A in the neutron star binary J0737-3039 will 
soon be measurable through detailed measurements of the periastron advance. We 
present the calculation of the moment of inertia of neutron stars with the masses of 
the components of the binary J0737-3039 for a broad range of equations of state of 
dense matter and discuss the implications of such measurement for constraining the 
equation of state. An observational determination of the moment of inertia of the 
pulsar A in J0737-3039 with the accuracy of 10% shall narrow down considerably the 
range of viable equations of state. We also show that limits on maximal mass of a 
neutron star provide a complementary set of constraints on the properties of dense 
nuclear matter. 
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1 INTRODUCTION 

Constraining the equation of state (EOS) of dense mat- 
ter in the interior of neutron stars is a fundamental prob- 
lem of modern astrophysics. At present, many different 
theoretical EOSs at densities exceeding nuclear density 
(po = 2.5 x 10 14 g cm -3 ) are available. Their predictions 
diverge at a few times po, which reflects poor knowledge of 
the EOS of neutron star cores. To constrain the EOS, one 
needs to measure neutron star masses and radii or some 
functions of these quantities. This is an extremely difficult 
task. 

Neutron star masses are best known for the radio pul- 
sars residing in binaries. The estimates of their masses have 
been achieved by measuring the post-Keplerian parame- 
ters of binary motion. The best studied cases are the dou- 
ble neutron star binaries, like e.g. the Hulse- Taylor pulsar, 
or the recently dis covered binary pulsars J0737-3039A,B 
feurgav et al J 12003:1. Masses have also been str ongly con- 
straine d in some binaries with white dwarfs (iNice et alJ 
l2004allbl) . Masses of the accreting neutron stars have been 
constrained by measuring the motion and spectral properties 
of the companions. The latter, however, suffer from larger 
systematic uncertainties because of difficulty in modeling at- 
mospheres of stars illuminated by X and gamma rays from 
the neutron star. 
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Constraining a second parameter like e.g. radius or 
any function of the mass and radius is much more diffi- 
cult. One can achieve that by measuring the thermal UV 
or soft X-ray radiation from the surface of a neutron star 
JPavlov fc Zavlirj|2003Tl . The estimates of the radius de- 
pends in this case on precise modeling of the spectrum in 
the presence of strong magnetic field and with unknown 
chemical composition of the atmosphere. The constraints 
on the radius depend also on the estimate of the distance to 
the source. Identification of spectral lines in the soft X-ray 
spectra would provide an estimate of the gravitational red- 
shift of a neutron star. Howev er, no such line have been re- 
liably identified yet (results of ICottam et al.ll2002l still need 
a confirmation). The discovery of kHz QPOs raised hopes 
to constrain masses and radii of neutron stars. These oscil- 
lations are most likely connected with motion in the inner 
disk around a neutron star. The estimates based on the QPO 
frequencies depend on a particular model of the QPO phe- 
nomena and also are uncertain due to erratic behavior of 
the oscillations. 

The discovery of the binary pulsars J0737-3039A,B pro- 
vided an opportunity to test General Relativity and pulsar 
physics in a unprecedented detail. The system is highly rel- 
ativistic, and only a few months of observations led to the 
constraint s similar to those of the Hulse- Taylor system after 
a decade llLvne et alJ 1200-3) . It is a nearly edge on binary 
which provides an opportunity to investigate the region in- 
side the light cone of pulsar B. The measurements of the 
pulsar motion are precise enough to measure the second or- 
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der post-Newtonian (2PN) effects. At this level of accuracy 
the periastron advance depends also on the spin-orbit cou- 
pling which involves the moment of inertia of pulsar A. The 
possibility of using this effect to mea sure a neutron star 
mome n t of inertia was firs t note d by iDamour fc Schaeferl 
(1988). lLattimer fc Schut J l)2004h estimated that the mo- 
ment of inertia of pulsar A should be measura ble with the 
accuracy of about 10%. iMorrison et alJ <l2004f) presented a 
calculation of the moment of inertia for three classes of equa- 
tions of state and showed that such measurement would al- 
low to distinguish between these classes. 

The moment of inertia of a neutron star is a very in- 
teresting quantity as it involves the square of the radius. 
Therefore, for a given mass, it is very sensitive to the de- 
tails of the EOS, and in particular its stiffness. Sensitivity 
of I(M) to the EOS was already utilized in the analyses of 
the energetics of the Crab Nebula expansion, powered by 
the Crab Pulsar at the expense of the neutron star rotation 
energy. Re sulting values of I were rather high, indicating a 
stiff EOS feeiger fc Haenselll20oi 120031) . This result is in 
accordance with very recent measurements of the mass of 
the PSR J0 751+1807, which a t the 95% confidence level is 
2.lt Q .t M llNice et al.ll2004al) . 

In the case of PSR J0737-3039 A the gravitational mass 
is already very precisely known (1.338 M©) and perspectives 
of the future measuring of I a via analysis of the effect of the 
spin-orbit coupling in the pulsar timing motivate a detailed 
calculation of the moment of inertia. Here we use a broad 
set of EOSs in order to study various correlations involving 
the type of the EOS and I a, and discussing observational 
constrains from measurements of stellar moment of inertia 
and mass. In section 2 we briefly describe the EOSs we use, 
and we present the method of calculating the moment of 
inertia. Section 3 contains the results, while in section 4 we 
present the summary and discussion. 



2 CALCULATING MOMENTS OF INERTIA 
2.1 Description of EOSs 

To illustrate the broad range of possible values of the 
moment of inertia, we use a set of twenty three EOSs 
of dense matter ava ilable in the literature. Similarly as 
IMorrison et al] <l2004f) . we divide the EOSs into three classes. 
The first class contains EOSs involving only nucleons and 
leptons, based on realistic nuclear interactions fitted to 
nuclear-physics data, and obtained using modern, precise 
solutions of the m any-body problem (Class I : APR of 
lAkmal et ai]ll998l W FF1-3 oflWiringa et al.lll988l. BB B1-2 
of iBaldo et al]|l997l. SLv of iDouchin fc Haenselll200ll . FPS 
of Pandharipa nde fc Ravenhalllll989l) . Out of these EOSs, 
only FPS and APR models were considered by Morrison et 
al. (2005). We added two EOSs based on the Brueckner- 
Bethe-Goldstone theory of nuclear matter (BBB1, BBB2), 
one EOS derived from the Skyrme effective nucleon-nucleon 
interaction, consistent with variational calculations for dense 
neutron matter (SLy), and three EOSs (WFF1-3) obtained 
in variational calculations of Wiringa et al. (1988). 

A second class contains EOS with less solid nuclear 
and many-body basis, but allowing for more speculative 
extrapolations : non-relativistic effective forces and rela- 



tivistic mean-field models involving nucleons and hyper- 
ons, and models with phase transitions to kaon conden- 
sates and quark ma tter (Class II: GN H1-4 of lGlendennind 
Il985l WGW1- 2C oflWeber et alJll99ll BGN1H2, BGN2H1 
BGN2H2 of iBalberg fc Gall ll997TlBalberg et alJ Il999l. 
GNHQ of iGlendennind Il997l (mixed-phas e had ron-quark 
EOS, model 2) and PREPL of iPons et ai]l2000l) . This set 
of EOSs is much larger and richer than Class II of Morri- 
son et al. (2005), who restricted themselves to two EOSs 
obtained using the Relativistic Mean Field model. 

Finally, third class is the most exotic, as it includes 
models of self-bound stra nge (quark) stars (Class II I: DB - 
DRS of iDev et alJll998l and SM, SM60 of IZdunikl Epoch . 
This class contains only three EOSs because the scaling re- 
lation (see Sect. 3) enables one to easily predict results for 
any EOS of strange quark matter intermediate between the 
stiffest SM60 EOS and the softest DBDRS EOS. 

Moreover, we considered additionally three specific 
EOSs of extremal c haracter. Two of the m are stiff (MF T17, 
lHaensel et ailll98lT) and soft (BPAL12, lBombacilll995l) ex- 
tremes of the nucleonic EOSs; they set two boundaries 
for the baryonic EOSs. While they reproduce experimen- 
tal saturation density, binding energy, and symmetry en- 
ergy of nuclear matter, they strongly overestimate (MFT17) 
or underestimate (BPAL12) nuclear matter incompressibil- 
ity at saturation. We also consider the CLES EOS that 
does not belong to any of above classes I-III. The causal- 
limit (CLES) EOS maximizes the maximum allowable stel- 
lar mass. This EOS assumes the SLy EOS below density 
2po = 5 x 10 14 g cm -3 and the stiffest possible (sound speed 
equal the speed of light) EOS at higher density. 

2.2 Numerical technique 

The moment of inertia of a neutron star is calculated as- 
suming rigid rotation, and in general case it is defined as 
J/fi, where J is the stellar angular momentum and Q is the 
angular frequency of rigid rotation as measured by a dis- 
tant observer. If fi is much smaller than the mass-shedding 
(Keplerian) limit for a given mass, f2 ms , then a slow ro- 
tation approximation is valid, and I can b e calculated us- 
ing non-rotating stellar model jHartlelll967l) . In the case of 
J0737-3039 A pulsar Q = 276.8s" 1 is less than one tenth 
of fi ms , and therefore the Hartle approximation should be 
valid within a percent (corrections to it are oc (fi/57 ma ) 2 ). 

The values of moments of inertia for hydrostatic stellar 
models based on the presented EOSs were obtain ed while 
solving the Tolman-Oppenheimer-V olkoff equation dTolmanl 
ll939HOppenh cnncr fc VolkofligH for spherically s ymmet- 
ric s tar, within the slow-rotation approximation jHartld 
1967). Additionally, for a representative subset of EOS we 
have directly checked the accuracy of this approximation by 
computing the moment of inertia for the axisymmetric rotat- 
ing stellar model within the highly accur ate 2-D LORENE 
framework feonazzola et alJll993L Il998l) . The virial crror 
indicators which describe the accuracy of computation were 
typically 10 -5 or less. As one could suspect, the difference 
between non-rotating and axisymmetric rotating models at 
this rotation rate is of the order of ~ 1% or less. It is 
principally dominated by the difference in methods of in- 
terpolating the tabulated EOS. Therefore, we draw the con- 
clusion that the spherically symmetric calculations are per- 
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Table 1. Highly accurate 2-D results for rotating axisymmctric 
stellar models of gravitational mass M = 1.338 Mq and rota- 
tional frequency Q = 276.79 s — 1 computed using the LORENE 
subroutines for a small subset of EOSs. The symbol p c ,i4 denotes 
the central density (in units of 10 14 gem -3 ), J47 [10 47 cm 2 gs _1 ] 
- the angular momentum, 745 = J/Q [10 45 gem 2 ] the moment of 
inertia and i? c i rc the circumferential radius (in km). In general, 
the values of the moment of inertia differ by approximately 1% 
from the non-rotating results - the difference is dominated by the 
error of the tabulated EOS interpolation method. 



EOS 


Pc,14 


J47 


I45 


-f^circ 


BPAL12 


19.57 


2.72 


0.98 


10.52 


APR 


9.74 


3.39 


1.23 


11.35 


SLy 


9.39 


3.57 


1.29 


11.76 


BGN2H1 


6.09 


4.53 


1.64 


13.47 


GNH3 


6.18 


4.77 


1.72 


14.28 



fectly suitable for reasonable estimation of the moment of 
inertia for J0737-3039 components. The results for the slow- 
rotation approximation for both pulsar components of the 
J0737-3039 system are displayed in Fig. and the precise 
parameters for a small subset of EOSs (for pulsar A) are 
collected in Tabled 



3 RESULTS 

The dependence of the predicted values of I a and 7b on 
the assumed EOS of dense matter is visualized in Fig. 1. 
Within the considered set of EOSs 7a, 45 = 0.7 — 1.9 and 
Ib,as = 0.7 — 1.7, a huge scatter indeed (we adopt the nota- 
tion 745 = 7/10 45 g cm 2 ). However, if we restrict ourselves to 
most recent EOSs consistent with experimental nuclear mat- 
ter parameters and involving only nucleons, then the range 
of predicted values narrows considerably, to 7a, 45 = 1.1— 1.3 
and 7b,45 = 1.0 — 1.2. The EOSs which involve hyperons 
and/or exotic phases of dense matter (quark and kaon con- 
densate) at supra-nuclear densities turn out to be so stiff 
for neutron stars with the masses of the pulsars A and B 
that they yield 7a, 45 = 1.4 — 1.7 and 7b, 45 = 1.3 — 1.6. 
Their high stiffness within PSR J0737-3039A,B results from 
the requirement that they should produce maximum allow- 
able mass > 1.44 Mq. The maximum mass is sensitive to 
the high-density segment of the EOS, with strong softening 
due to hyperons or a phase transition. For this segment to 
be sufficiently stiff to respect consistency with observations, 
the lower-density segment (without hyperons and/or new 
dense phase) should be indeed very stiff. The two extremes 
of nucleon EOSs set are MFT17 (very stiff) and BPAL12 
(very soft); while they reproduce experimental saturation 
density, binding energy, and symmetry energy of nuclear 
matter, they strongly overestimate (MFT17) or underesti- 
mate (BPAL12) nuclear matter incompressibility. 

The values obtained for the causal-limit CLES EOS 
(open star symbols) are similar to those obtained for the 
Class II EOSs. At relatively low neutron star masses, the 
SLy envelope with p < 2po is important for stars structure. 
However, the core with v S ound = c shifts up the values of 




r ✓ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1 1.2 1.4 1.6 1.8 2 2.2 



M [MJ 

Figure 2. The dependence of the moment of inertia on the grav- 
itational mass for a few representative equations of state. 

7a, b, by more than twenty percent, compared to the SLy 
ones. 

A separate subset of EOSs (open triangles) is composed 
of three models of self-bound strange quark matter. The val- 
ues of 7a, 45 and 7b,45 predicted by the strange-star models 
are very scattered 0.7—1.7. However, this range of values can 
be easily understood in terms of the scaling properties of I 
and M with the value of the surface density of strange stars. 
To a very good approximation the values of 7a and 7b calcu- 
lated for different EOSs of strange quark matter are related 
by a simple relation involving surface density of bare strange 
stars (= density at zero pressure), p s - Namely, 7a calculated 
for the EOS with zero pressure density p s is related to I' A> 
calculated for the EOS with p' s , by I' A ~ (p B / p' s ) 2/3 Ia (this 
relation is obtained using constant-density Newtonian star 
models). Therefore, huge disparity in 7a for the DBDRS 
and SM60 EOS can be accounted for by the difference in 
p a - We considered bare strange-star models of J0737-3039 
A, B. However, the mass of the crust of strange stars is at 
most 10 -4 Mq, and therefore its contribution to 7 is at most 
0.01% and can be safely neglected. 



4 SUMMARY AND DISCUSSION 

We have presented the moments of inertia of neutron stars 
J0737-3039 A, B for 25 different EOSs of dense matter. In- 
spection of Fig. 1 leads to conclusion that a measurement 
of 7a with a 10% precision would necessarily imply rather 
precise picking up of the EOS of matter in J0737-3039 A. 
It could be a stringent test of "realistic EOSs" of Class 
I. Either the EOS inside J0737-3039 A is stiffer than our 
best nuclear models predict and belongs to Class II, or is 
in harmony with nuclear many-body theory represented at 
its bes t by Class I. We there fore confirm and extend the re- 
sult of lMorrison et all i2004l) . While it would be difficult to 
eliminate strange star model altogether, measurement of 7a 
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Figure 1. The moments of inertia of the binary pulsar components: A- left panel, and B - right panel. Filled circles: Class I (realistic 
EOSs involving nucleons and leptons). Open circles: Class II (relativistic and non-relativistic mean- field models involving hyperons, kaon 
condensates and quark matter). Open triangles: Class III - strange quark stars; the range of values of the moment of inertia available for 
strange stars lies between the triangles labeled DBDRS and SM60. Open squares: extreme nucleonic EOS. Open star: causal-limit EOS 
(CLES). For further explanations see the text. 
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Figure 3. Moment of inertia of a neutron star with mass of pulsar A versus the maximum mass of a neutron star plotted for the same 
set of EOSs as in Fig. 1. The vertical thin solid line denote the 68% confidence lower limit on the maximum mass (2.0 Mq), whereas 
the thick solid line marks the 95% confidence limit (1.67 Mq) - the limits were obtained from the PSR J0751+1807 mass measurements 
(see text for details). 



would strongly limit the value of p s of strange quark mat- 
ter at zero pressure. The measurement of the moment of 
inertia 7a will be sensitive to the properties of the EOS for 
the densities typically in the range of the central density of 
the pulsar A. From Table 1 we see that this corresponds 
p < (6 — 10) x 10 14 gcm~ 3 , i.e. or p < (2 — 4) x po where po 
is the nuclear density (see Sect. 1). In this estimate of the 
density we have neglected the unrealistically soft BPAL12 
model. 

It should be stressed that hopefully stringent con- 
straints on Ia do not probe the EOS at densities crucial 
for the value of maximum allowable mass for neutron stars, 
-M max . A good example of two EOSs which yield very simi- 



lar values for I a and Jb but very different stars at M max are 
BPAL12 and WFF1 (Fig. 2). However, as Fig. 3 shows, ex - 
istence of 2.1 Mq < M max neutron star iNice et alJl2004al) 
combined with determination, within 10%, of Ia, would 
quite precisely establish the complete EOS of dense mat- 
ter inside neutron stars. Lower limits on the maximum mass 
of a neutron star are sensitive to the properties of the EOS 
at the density p « (6 — 10) x po, much larger than densities 
relevant for the constraints from measuring I a- Both con- 
straints together would constitute an ultimate test of the 
many-body theories of dense matter. 

While the measurement of the moment of inertia of the 
pulsar A in the J0737 system is approaching it is interest- 
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ing to speculate about further investigations of the neutron 
star equation of state. Measuring the moment of inertia of 
the pulsar B in this system does not look feasible with the 
current instruments. However precise monitoring with the 
next generation instruments like e.g. the Square Kilometer 
Array may yield highly accurate data that enable determi- 
nation of 7b , or moments of inertia of neutron stars in other 
binary pulsar systems. The Parkes Multibeam Survey may 
also bring discoveries of other highly relativistic binary pul- 
sars, and provide measurements of the moment of inertia 
of neutron stars with other masses. One should note re- 
cent measurements of masses in neutron star white dwarf 
binaries which hint at the existence of neut ron stars with 
masses as large as 2.1 Mq jNice et al.ll2004 a|[b|). Such mea- 
surements provide lower limits on the maximum mass of a 
neutron star and thus are complementary to the measure- 
ments of the moment of inertia. Limits on the neutron star 
maximum mass probe the densest segment of the EOS. We 
illustrate the strength of the constraints coming from the 
measurement of the moment of inertia of pulsar A in J0737- 
3039 system together with limits on the maximum neutron 
star mass in Fig. [3] The two values - I a and M max are not 
correlated since they probe different ranges of densities of 
the EOS. A simultaneous measurement of the moment of 
inertia of pulsar A with accuracy of ~ 10% with a limit 
on the maximum mass of neutron star approaching 2.0 Mq 
should narrow down the possible set of nuclear matter EOS 
to just a few possibilities in the near future. We also plot 
the lower limits on the maximum mass of a neutron star 
implied by the recent measurement of t he mass of the pul- 
sar PSR J0751+1807 jNice et al.ll2004al) . We have obtained 
the limits shown in the Figure by assuming that the errors 
on the mass are approximately gaussian. We have then in- 
tegrated the approximate probability and obtained the 68% 
lower limit on the maximum mass of 2.0 Mq, and the 95% 
limit of 1.67 Mq. These values are plotted as vertical lines 
in Fig. |21 Thus a simultaneous measurement of the moment 
of inertia of pulsar A with accuracy of ~ 10% with a limit 
on the maximum mass of neutron star approaching 2.0 Mq 
shall narrow down the possible set of nuclear matter EOS 
to just a few possibilities in the near future. 
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